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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE NO. 1002 


STRESS ANALYSIS OF COLUMNS AND BEAM COLUMNS 


BY THE PHOTOELASTIC ' METHOD 
By B. F, Ruffner 

SUMMARY 


Principles of similarity and other factors In the 
design of models for photoelastic testing are discussed. 

Some approximate theoretical equations, useful in the anal- 
ysis of results obtained from photoelastic tests, are de- 
rived. Examples of "the use of photoelastic techniques and 
the analyses of results as applied to uniform and tapered 
beam columns, circular rings, and statically indet erainat e 
frames, are given. It Is concluded that this method is an 
effective tool for the analysis of structures in which 
column action is present, particularly in tapered beam 
columns, and in’ statically indeterminate structures in which 
the distribution of loads in the structure's Is influenced by 
bending moments due to axial loads in one or more members. 

• ; ; • •!-. 'i i 

'Introduction • ' '■ 


Many authors have discussed methods for determining the 
bending moments and buckling loads in columns and beam col- 
umns. Niles and Newell (reference l) give solutions of the 
problem for various loading conditions v/ith the El of the 
beam constant for the span under, co_ns.ide.rat ion . Timoshenko 
(reference 2, p. 128) gives the critical buckling loads for 
tapered columns of several types. . The problems discussed 
in those references represent solutions of the basic differ- 
ential equation of a beam loaded with an axial load P and 
lateral loads q (x) . This differential equation is: 


. ;r 

VT A ; 5 

2 j x * 

dLx 3 


+ Py = f(x) 
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When 31 is a constant, this equation- yields readily to 
integration and the bending moments at any section x 
are given by: 


M = *f (x) - Py ( 2 ) 


The critical load P is: 

cr 


cr 


Ctt'EI 

l 2 


(3) 


where C 'Is the fixity coefficient- depending on the end con- 
ditions of the beam. The value, of . P is independent of 
the. lateral loading q(x). .... . cr 

In many cas : es, however, El is not constant “but varies 
either continuously or. d is.c ontipuo.u^ly over the length of the 
span under consider r^ti on., Equat'd on ( 1). may then be integrated 

"by approximate numerical metho.ds/ .General solutions are 
possible only in a .£ e»w c.ases. 

• * 

It is the purpos.e of this . ijqve s t i gati an to develop a 
simple,- accurate, and quick method, -of determining moments 
and critical Loads in members loaded , as beam, columns by 
the use of methods of photoelastic analysis. The method 
will be most useful for statically indeterminate structures 
containing members with nonlinear behavior. 

This investigation, conducted at the Oregon State 
College, was sponsored by and conducted with the financial 
assistance of the National advisory ComE'.ittee for Aeronautics. 


DESIGN OF PHOTOELASTIC MODELS 
General Considerations 


In order for a model test to be useful it must be 
possible to predict accurately results on the full-scale 
structure from results of the model test. To accomplish 
this purpose, it is necessary to design photoelastic 
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models with several fundamental. principle. s in mind. First, 
principles of similarity must* be met between the model and 
full-scale. Sec ondly, : .• t he mode 1 ,must be constructed so that 
limitations imposed by techniques of testing are considered. 
Thirdly, the models must be of such dimensions that they can 
be constructed within -.tolerances .required, 

* - -t • r . • <7 : • r . r j. .. t 

Similarity Principles:' 

By: use of rlimen s i onel nna lys i s r i t is possible to de- 
termine the conditions of similarity that must be met in 
structural models. (See reference 3.) If P is any one 
external load on a structure, , V ;1 P,V. fc P V 3 P. . f , are any 
other external loads on the structure. i. . li , t, . I 4 . . . . 

are pertinent linear dimensions of the ‘ structure, E is' 
modulus of elasticity, p is Poisson’s ratio, then it may 
be shown that, in general, 


R 


= f 


.(A 

\3l ^ 


> V 1 »*V 2 **V 3 


l z l 3 l 4 

» ” > » ” 
1 i 1 2 . ^ 1 


( 4 ) 


- 1 " id i V - i ■■■ 

;J .1 


* ; .M . f P; ■*'> i - .1 ..,l a l 3 - 1 4 - A 

. - f 2 ( ^ i 2 ‘jY 3 I * t. . i. -T— , T-— -r~ ... |4 . J ( 

•Pii • ■ \ . I J x ^ i a ; / 


5) 


( 


where R is any reaction or ,,i n.t ern.al force, ’in a; member, and 
M is a bending moment in a member. * •' 


If all the parameters on the right-hand side of the 
equatioh are the same’ on the” model ahd the -/full— scad e struc- 
ture, the left-hand members will be the same for the model 
and the full-scale structure. Values for the full-scale 
structure would then be obtainable from model test. 


In certain particular cases it.,is, unnecessary to have 
all the nond imenai onal parameters the samb /on the model 
and the full-scale structure, since, some are not pertinent. 
Several types of structures will be briefly discussed. 

: ’ • ’ •’ *•' * '* d V : it \ . •. ; . * t * * . 
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i- STRUCTURES IN ' WII CH 'MOMEHT-S , ' INTERNAL LOADS, REACTIONS, 
AND SO FORTH, ARE LINEAR FUNCTIONS OF THE 'LOADS 


In many -structures the forces in the members, the 
stresses, deformations, and so forth, are all linear func- 
tions of the load. Then the nond imens i onal parameter 
P/El, r is not significant. If "the ratios R/P and M/P^ are 
determined for any one set of loads, they are then the same 
for any other set of loads in constant 1 ratios to the first 
set. ' • ' 

In .most structures' of. this type Poisson’s ratio |* is 
not significant either. 'For instance, consider the prob- 
lem of .bending of beams in whi-ch-’the ’equation: 


1 2 

d* 

is applicable. Since p does not appear in this equation, 
it may be considered permissible to neglect this parameter. 
If the beam is loaded so that M is not a function of y, 
then both P/El^f and p may be neglected. Furthermore, 
geometrical similarity is. not ne cess pry for the determin- 
ation of moment s t reactions, and so forth , since the linear 
dimensions of the cross sections come in only as influencing 
the moment of inertia. Consequently, the similarity con- 
ditions- to b-e met wuv.il d "'be' 'expressed' by the equat-ioiis : 




l 3. 

1 






,v 2 ,V 3 



1 1 1 1 




( 7 ) 


These conditions require that the loads on the model and 
the full-scale structure be in the same ratio, that the 
lengths l lf 1., I3 . . . of the neutral axes of various mem- 
bers of the model and the structure be in the same ratio, 
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and that the moments of inertia of various parts of the 
structure he in the same ratio on the model and the full- 
scale structure. 


COLUMNS AND B DA M COLUMNS 


WITH STATICALLY DETERMINATE ■ 


EXTERNAL LOADS 


Structures which behave nonlinear ly, hut which are 
loaded so that the material is in the elastic range every- 
where, include those structures classed as columns and beam 
columns . 

The basic equation of the elastic curve for these is: 


El 


dx 2 


M 


but in this case M is a function of both x and y. For 
instance, consider a beam column of length l with 31 = 
f (x) . The equation of the elastic curve of this may be 
written: 

' \ V ' + pi V y ^ 

" d Y x\ s SI M ^ 

where P is the axial compress! 

the bending moment due. to lateral loading W, composed of 
loads y iP , y 2 P . and . so ■ f or t'h-. ' ’ 

In order for this equation to be identical for both 
model and full-scale structure, the conditions that 


Wl- „ 

= k 1 i . . 

SI V l 



ve load and kW f M — ) is 
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X 

must "bo-met at every point defined by - . The subscript m 

l 

is used to denote model, and the subscript s to denote 
full-scale structure. Dividing the first of these equations 
by the second gives: 



Here ajgain Poisson's ratio p, is not pertinent. The bending 
moment may then be written: 


M 

Wl j 



I 2 l 3 P 
» * — / 
l a l 1 31 / 


( 8 ) 


It is usually more convenient to let the moment, at some 
point, due to all lateral loads be M 0 and write the fore- 
going equation in the ai te mate fora, 


^0 

M 



^ 2 ^ 3 P^l 3 ^ 

• • * . » I j 

l 1 31 / 


(8a ) 


If I is a variable, this must be written: 



(8b) 


where I Q is the moment of inertia of some selected section. 
The conditions for similarity for a column or ho am— column 
test are then : 

1. The ratios between loads on the model and between 
corresponding loads on the full-scale structure should be 
the same. 
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2 . The ratios of I Q /l f or qpTresponding' p'oints in 
model and full-scale structure should be the same, 

3, Ratios of lengths of neutral axes on model and full- 

scale structure should.- be. the same; : l.\ 


4. Values of K 0 /M then will be the same model and 
f ull • s ca^e-. if the patio P \ ?/BI 0 - is t he >same and full 

s c a 1 e . ., 1 , . . j . , ... , % ■ ( f - i ■ j • * : • ■ 





•; o 



STATICALLY I DETERMINATE . STRUCTURE $.- IN WHICH ONE-, OR MORE 



MEMBERS 'IS A COLUMN OR A BEAM COLUMN 


l 7 ’ 


O i 


a ? . ’ ? [■ 


-i ! 



- ;* Ci l '■ -V: ■ ■ 1 t' ... : ‘ ' Xt ' '» ’• * O is Si ’> ”• •*' 1 

. ; For. structures of this type all similarity- c^ndfitions 

applies, hie to., teau columns must be met . addition, ..‘h.ow— 

ever, the, d i st r.i but i qn qf loads may depend, on the cross- 
sectional areas of trie 'members as well as tHe.ir moments. ,. of 
inertia. Ih structures in which it' is known, t,hat 'de.fo.rmr- 
atioris due to direct t’e‘nsil l 6 and compressive .st r e's s e s. „ar e 
unimportant , it will' be .sufficient t c :..eot on Ly . thos e simi- 
lar it y conditions listed, for .beam columns, If', however,' 
both bending' 'dVrbrnat ion's. ' arid d. ef orma.t,i o.ris due. to tens il’e... 
and: compres si ve stresses ar e impor tant theri t.he cpnd it'i on, 
that the areas of all members must be in the same ratio 


model and full scale is important. Equation (5) then may 
be written in the fern:' 

, •’ ,» ! . :.S 


t t\ ' 





I ■ ' ‘k e/ 


, -° 


U 

A Vr 




A 0 /A is tiie ratio of the area of some selected section of a 
s t rue tur al r member t o thd area of any other, section^ If; ' ( 

structures • : of the type in’which ’shear de format ionb "p l'$ p\ 

are important, then the conditions 'of ge ome t r i c ! s im v i iar i t y v ‘ 
•must be met in addition to the foregoing. This in most 
structures, however, is net the case, and the model need not 
be geemet r ically similar s te full scale. Usually models with 
members of constant thickn.ess may be* used to represent full- 
seal e c ond i t ions* 


8 


NACA TN No. 1002 


DESIGN OF MODELS IN WHICH ONE OE MORE MEMBERS ARE LOADED 
-AS BEAM COLUMNS BUT IN WHICH LOADS IN MEMBERS 


If this structure is of the type in which deformations 
due to tensile and compressive forces in the members, are 
negligible as compared to deformations due to bending, then 
the method of design of the model and the determination of 
test load conditions are the same as for the beam-column 
models previously discussed. If, however, deformations due 
to tensile and compressive loads may be important, then it 
is necessary to meet the additional similarity condition 
that areas of various members on the model must be in the 
6ame ratios as corresponding areas full scale. In addition, 
the ratio P/Ea q for the model and the full-scale structure 
must be the , same . To meet these conditions, in addition to 
the other conditions of similarity, necessitates a variation 
in thickness of! the model material. If a structure is com- 
posed of several members each of constant El and constant 
area, the model may be constructed with little .difficulty if 
the thickness. of the phot oela st-i c model is different for 
each, member...' The model -design and loading conditions then 
will be determined.' as- f o>tl ow s'? 

1. Assume some convenient' over-all dimension for the 
model. Let this be L om . 

2. Next designate the moment of inertia of some full- 
scale member as I Qg and the area of. that member A QS . 

Choose a moment of inertia which will give a reasonable model 
dimension, say I Qm , for the corresponding member on the 

model.. The dimensions of this section of the model member 
then may be determined to meet the necessary similarity con- 
ditions for area. From equation ('9 )•» ... 


ARE STATICALLY INDETERMINATE 




and 
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Dividing the first of these by the second gives: 


ftoVS „ 
v *• 


A 0 1 » 


f“) 


or 


a_ = a _ _ Ids 


■om 


os 




( 10 ) 


o s 


If I Qm has been chosen, then the dimensions of a rectan- 
gular section of the model may he determined. If b is 
the thickness of the model and h the v/idth of the section, 
then, 


A = b^h„ 
om oo 


om 


»o h o 3 

12 


Substituting this in equation (10) gives: 


b o h o “ A os 


L --\ 2 b„hj* 


OS' 0 0 


I 


om 


•' 121 


0 s 


- ^ ~om\ 

v ' L I 
o s 


121 


0 s 


0 s 


> 


J 


(11) 


Then, 


■ (r 1 ) (r 2 ), 

om os 


‘os 

121 


0 s 


( 12 ) 


10 


NACA TN No. 1002 


and the model dimensions are determined for that section. 


To determine other model dimensions, say at a section 
denoted by the subscript 1 , the conditions, 


^b 1 h 1 ^ 

Vb 0 h 0 V 



must be satisfied. 



-Vb' h 

xi ■ ) . 

u o n 0 

■ • , • . } 

0 s 

These give, 




(11a) 


(12a) 


In some cases it may be necessary to try several dif- 
ferent values of I before a, model of convenient dimen- 
sions is obtained. For the model representing structures 
composed of members each of constant area and moment of 
inertia, the model will then have its members each of con- 
stant thickness, but the thicknesses will be different for 
different members. 

V/hen structures have individual members of varying 
moment of inertia and area, both the width and the thick- 
ness of the members of the photo elastic model should vary 
in accordance with equations (11a) and (12a). It is not 
practical, however, to construct photoelastic models in 
this manner. Usually, the condition that the moments cf 
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:.v.z rs'.'.s 'of tii 


he .models must vary according to f—^) = (—') 

■ : cl I. o v : vrf \I 0 / VI 0 J 

■ ■ " • • ’ .% - • m - . E 

i nri or ~ k n f o n nti tr a i i a 1 1 Tr ... -m .n tr "h . /r> "h _ 


-s the most important. Good accuracy ilsually-.rfflay beobr- 
tained if the thiclrness of a member is determined so tha.£, ff _ 

^ A av ^ _ /A 


its average area will meet the. conditi on 


- 

•o ; m VA o 


m w s 

The correct variation in moment of Inertia is then obtained 
by varying the width of that member . 

3. When the dimensions of thermodel have been deter- 

, ‘ \ • ' •* ' ~ , 

mined, the loads to be a^Pl ied t 0 ; the. model tijqn are deter- 
mined ’ from the condition that'- 


a > 
’ ? 


- ..•? •> i - :■ 

. , , m t * c 


/.fh 0 2 h; ; = 00 s 
v S I 


’ J 

t r 


l 


31, 


o? A :■ 4 

) * 0h.< ,- s . 

BT / 


m 

*< , : .i ; : 


i f ~ tch 


ADDiTEGSrAL ' GOMMSHT.S pll MODEL DESIGN 


o.- f , 




. ,l,.,The model should be large enough so -that reasonable 
t- olei ance's in ths- cpms-tru-chipp pf.the model , will not cause 
seripus errors i n ’ the.- v.a.ljie s .,of . the momeht of inertia.. ?o 


instance, 


pS.tjp.fi moment of inertia,.- for 
if the- mode 1- ca-n. be. cut to dimensi on 1 wi.t hi n ■ £ 0 . 002 

err o r 


in m’cine'nt' of 


model i's Cl. 100 ‘inch', 
give' sin error 


inertia qf 
If buck- 
of.' 6 percent 


•1 

5 r 


ir.cn, there 1 will be a possible 
±6 percent if the width of the 
ling loads are desired, this would 
in e st i mat ed -huciC'l ing >l ; qra-d. ,. r ■ .. 

j ■ • i • ■ i ■ $ • • 

i'. 'Usually, bdi t er fringe patterns, less 'f edge .'effect , 11 
and tetter photographs will be obtained 'with Bakelite of 
0 .300- inch* thickness • or less.,. If, the thickness is ’too small, 
however ^ t oo few fringes may be obt ei ned , f or the r all owabl e 
stress in the Bakelite. The accuracy of the tests., particu- 
larly i ; h' : ' determining., bending, moment s and axial loads, is - im- 
proved if the model is loaded so as tp produce a large number 
.,of fringes. Thicknesses of Bakelite BT 61- 89’3 j T r on* 0.. 125- r 
inch .to' 0.3'00 iweh& .give .good- results. , 

• \ ••• •' A ' r ! , ,, ■■ ... . , ‘ ' ' ! * X -> r < ; .... 

.... • r • '-'h'-m .« ri .) >■ . 

' , . '3 , !, Tn '"Cdhst r-ticti ! hg model s, ,' ( of statically l'nd:eterminat e. 
st J T‘u ; c'’tur'e s : i t' -is well t ^ o -ke-ep , i,n mi nd that slight /ina-ecara--. 
cies in lining up supp or t s , . f it' t.ing member s ' t d'geth'er , -a-nd. .so 
forth, will cause internal loads ‘in thh structure swh-en +no ^ex- 
ternal loads are acting. At times it is desirable to deter- 
mine the effect of misalinement of the supports on the stresses 
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This, of course, may be done by following similarity princi- 
ples outlined previously in determining the correct model 
mi salinement . 


PRELIMINARY TESTS 
Purp o se 


.• : The purpose of these preliminary test6 was to inves- 
tigate experimental procedures for determining bending 
moments and critical loads in pin-ended columns and beam 
columns. Models were chosen which could be analyzed by the 
methods given in references 1 and 2 so that a measure of 
the accuracy of the experimental results could be obtained. 


Models and Tests 

Models were milled from sheets of polished and annealed 
Bakelite BT 61-893. Dimensions of Models 1 to 11 are given 

in figures 1 and 2. 

% • , 

• • Si. 

Model's "were '.Loaded in the frame shown in the photograph 
of figure ’3 i>-In -figure 4 is shown the arrangement of the 
loads " l o'n the models. The loading frame was designed so that 
it- 'wsis -possible t'O : 'ya,ry -the end load P . and the lateral 

load C4 independently* w - ; . . * . '*• - ' 

* > ,• • 

- ; . 

For each loading condition, a fringe photograph under 
circularly pb-lhrized J-ight vae takep-, of the center portion 
of the beam . A monocromatic light p f 5461 angst rom uni t s 
was’ used' in- All? test s. f Fringe photographs were numbered by 
the f oil bwi ng; system ; Phot o No. - Model No. - Test No. - 
P in l : b Q, -In lb. Figure 5 shows a series of fringe 
phot ograph-s ‘obtained on Model 6 with Q = 4.18 lb and P 
Varying from ; 0 .'O' to 30.0 lb. In table I tests run on Models 
1 to j il' ar e- '-LL st ed . • 7c . • j 

• ^ J 4. P "I ; , - 

Model 12, shown in figure 6, was used as a tension te6t 
specimen- '-for thu det er'mdn'ation ; of t h.e . mp4svl,p sl „of elasticity. 
This specimen was loaded up. to 3Q00 psji., The strains were 
measured by Eugg'enberger t ensomet er s . The modulus of elas- 
ticity was found t:0' be 668,000 pel'* 
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fie suits' 

;V *• .* ' ! ’* • ' ' ’ * * 

•• Because of the symmetrical loading on the models of 
tlfis preliminary; -series of test s, the sections at the cen- 
ters of the spans of the beams are all sections of maximum 
moment. At these sections of maximum moment the vertical 
shear is zero and the only stress, acting 'is the fiber stress 


a 


^ where 

I A 


is the distance from the neutral 


axis of the beam, measured in the plane of the loads, posi- 
tive toward the upper edge of the beam, and M is the 

c 

bending moment at the center of the span. Consequently, 
the Mohr 1 s circle of stress for any element in this cross 
section is as shown in figure 7. The principal stresses on 

this cross section are th.en <j. - - —9.1 ^ £ and a = 0. 

1 I A 2 

Since the fringes determined phot o elastically are loci of 
points at which the principal stress differences are a con- 
stant, the CTj may be determined immediately* from the 
fringe patterns if the fringe orders are known. 


The stress- op tic law may be written in the form: 



(13) 


where cr 1# cr 2 < are principal stresses, n fringe order, k 
a constant depending on the wavelength ,of. light, the prop- 
erties of the material used, and .t he., ^thickne s s of the model 
in the direction of propogation of the light. At 


z = - 


y. 

let n = n,p, and at’ z = - - let •• ji = 


n-j, . Then, 


X T 


v V " -M fe h * : P 

= knrp .= -r — S_~ r 

. * . V, 21 ■ .a .... 

1 • !. V. V. 


M h P 
ct i-r = . kn B = *- J — 
^ • •‘ 21 ’ A 


Solving these simultaneously for M gives, 

c ^ 


Ik 


M c « - (» T - n B ) 

h 


. .t i ■ ■ 


(14 ) 
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If the loading of the “beam is such that th r e bending 
moment is known from principles of statics, equation ( 1 4 ) 
may be used to determine k. In many tests the ratio M oc /M c 

is the desired variable, . From equation (14) this is: 


M oc _ ^ n T ~ n 3 ' P--Q 

- n 3^p 


(15) 


In table I are recorded the. fringe order differences 
(ns - n<r ) for all pr el iminary ,.t e st s run. Also in table I 
are given values of These were obtained by divid- 
ing (n .3 - foir .the test with no axial load; by the 

values of ( ng - nj') obtained on each of, the other tests of 

the same model' yith, varying values of the axial ■ load P. 

For Models 1 ..t. o 6 the theoretical critical loads are 

P cr = tt s 3l/l 2 . Tlpese were computed and entered in table I. 

For the taoered Models 7 to 11, the critical loads were com- 
puted by : i nt era olat i : ng the -table &: given by 1 Timoshenko in 
reference 2, page 128. These also are listed in table I. 


In figures 8, 9, and' 10 are plotted the results of 
tests on Models 1 to 11. Faired curves were drawn through 
the experimentally determined points and terminated at the 
computed values of P. cr . . -an inspection of these curves 
shows that points do not’ lie on th,e faired curves within the 
limits of error in the determination of moments' fr'oa the 
photoelastic fringe patterns. The scattering of these points 
was believed to be due t o, -?r i ct i.on in the end ' support s . 


An examination of these curves also shows considerable 
variation from the straight-line plots that were expected 
from theoretical considerations. The curvature of these 
plots was believed to be due to eccentricities in the appli- 
cation of the end loads; To check this, the following ap- 
proximate analysis of the effect of eccentricities was made. 


Approximate Theoretical .Analysi s of Effect 

of Eccentricities and End Couples 

If eccentricities existed on Models 1 to 11, it was 
known that these were small. The exact amount of these 
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eccentricities could not be determined accurately by measure- 
ment. Therefore, the following approximate analysis was made 
to attempt to evaluate these eccentricities from the test 
results. 

Timoshenko (reference 2, p. 82) has pointed out that a 
Fourier series may be used to represent accurately the de- 
flection curve of a beam column. He also shows that, if 
only the first term, of the series is used, the deflection 
curve, for many types of beams, is represented within re- 
auired engineering accuracy. In figure 11 the notation is 
defined. Then, let 


y = -a, sin — ( 1 6 ) 

1 l , . , 

■ • • W ; - ; . • , ; • ■ ? ■* 

where ~& 1 is the deflection at a point midway between the 

points of inflection. ‘The bending moment M, at any point 
x is, ' 

M = HI iiz - a, sin 22 (17) 

dx 2 r i 


In general, El may he a function of x. Let ("El ) ^ 

be a constant which may be substituted in (l7)'"to give the °' 
same value of the critical load that will be obtained if 
the actual variati on of SI is used. Equation (17) then 
become s : ‘ '• 


M 


n s (EI ) eff 


TT X 

sin — 

l 


(17a) 


The term 
column . 


(®*) e ff .will not be the average El 
It is a value which will satisfy the equation, ^ „ 


for the beam 

0 ’i' 




tt 2 (E I), e f f - 


cr 


It may be considered, to be the stiffness of a hypothetical 
uniform beam column having the same buckling v lo ad as the 
actual beam column under con si d era t i on . 
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Referring to figure 11, write the bending moment at any 
point xs 


M = M Q + f M 2 

L 


- P(y x + 6 x ) 


1 /L - 


a - x 


L • 




+ - P (y 2 + 6 2 ) 


la 


- Py 

L / 


(18) 


Consider conditions at the point midway between the points 
of inflection, at x/ l = 1/2. If the bending moment at this 
point is designated as M c , and the nortion of the bending 
moment due only to the lateral loads q(x) as M oc , then 
equating (17a) and (18) gives: 


M + 
oc 


ai 


Mi 


- * a,)] + 


Tt 2 (SI ) ef f 


(19) 


ana 




- P(y x + 8 X ) 




L - a l 


) 


“■ - p(y * + •* > ]6 ; &)} 


where 


p _ V(BI) ef f 

l 2 


( 20 ) 


( 21 ) 


P or certain loading conditions, equation ( 20 ) may be 
great ly . si mpli fi ed . Consider, for instance, the following 
example: Beam column with pin ends, no end couples. For 

this case, l = L, and = M s = 6 i = 6 3 = a = b = 0 . 

Then equation ( 20 ) becomes: 


M c = 


0 - P 


M - P 
“oc 


(^) 


or 
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0 - P p P 


r q lizi ±_JLal1 

P 

[l - l el 

2M 0C 

1 

8 

L . . . , H J 


where ; ■ ’ * ’ 

e = P(y i + Isl 

• " i 2M 0C 


( 22 ) 


(23) 


Let a = P/0; then equation (22) becomes. 



M c 1 - <xe 


( 24) 


•Squation (24) indicates several int er e.st ing,. re sul t s : 

r •* *V / . r ' ' ' 

1. If the eccentricities y x anc^ yg , are zero, then 

e = 0, and ^ oc /^ c versus P should plot as a straight 
line . 


2. The greater the value of the ( mo.me.nt. M oc due to the 
lateral loading', the less the ratio, j ? a,nd the 'less will be 
the- effect of any eccentricity on .-t.her rje sul't.isV' " " 


oi' 


3. If the plot of M oc /M c versus '+ P is' hot a straight 
line, ho'wevpr’j the experimental cur v e ’ fiu\y~ b‘e analyzed to de- 
termine the’ amount of eccentricity preseiit during the test. 


If the average eccentricity ZjL — ~ is 

2 : ;;; , „ 

equation (?22) may be written in the form: 


den ot ed 


Ly 



, M oc - M 0 v + Me 
P • 0- 



y av * 


( 25) 

i - 


The test results give values of M oc and M c for various 
loads P. Consequently, the experimental results may be 
substituted in equation (25) and (3 and y av computed. 
Por instance, 
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3*0C> 

and 

P x 

be 

the 

c or re sponding 

axial 

1 oad 

0 c » 

and 


be 

the 

corresponding 

axial 

load 


Substituting these in equation (25) and solving the result- 
ing equations simultaneously, gives: 


P 


1^2 
- p 2 


( 26) 


and 


y av 


!oc + 

*1 e 


( 27) 


Critical loads for Models 1 to 11 were computed from 
equation (26). Values of Pi and P 2 were taken from the 

faired curves of figures 8, 9, and 10. In table II are 
shown results of these computations. With the exception of 
tests on Models 1 and 2 the values of critical loads obtained 
from eauation (26) are in good agreement v/ith the theoretical 
values . 

The calculated values of the eccentricities, obtained 
from equation (27), indicate that relatively small errors in 
model construction will account for rather large variations 
from straight-line plots shown in figures 8, 9, and 10. 

Prom this analysis of the results obtained in the preliminary 
tests, several tentative conclusions were drawn in regard to 
the best test procedure. 

1, The factor e = ~ should be as small 

M 

oc 

as possible to avoid effects of eccentricities. This may be 
done by using rather large lateral loads. which make M QC 

large, or by keeping the eccentricity y av as small as 
possible . 

2. Since any error in M will appear in all points, 

O G 
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the fringe values for the leading conditions with no axial 
load should he determined to. a high degree of accuracy. 

Less possibility for errors will exist in the determination 
of M oc if the fringe order difference (nij - ng ) for this 

condition is in the order of 8 to 10 fringes. 

3. The model should "be loaded with axial loads up to an 
amount producing M c = 4M QC (approximately). This will 

assure a more accurate determination of values P^ and P 2 
for use in equation (26)., and consequently a more accurate 
determination of B. 

4. At a moment M q = 4M oc the stress a = -P/A - M c h/2I 
in the model should not exceed the elastic limit of the model 
mat er ial . 

5. The pin-end connections should he in good hearings to 
prevent 1 the effect of the formation of a friction couple at 
the pins. Pin friction will have the effect of givi.-ng- an 
effective increase in the fixity coefficient. 

TESTS TO CHECK CONCLUSIONS BASED ON RESULTS 
• - ' OR PRELIMINARY EXPERIMENTS 

Beam Columns with Eccentricities 

' • • v • .,tj' . 

Models and Tests 

: Models 13, 14, 15, and 16 were constructed as ‘shown in 
figures; 1, 12, 13, and 14. During tests the effect of fric- 
tion in the pin ends was eliminated as much as possible hy 
vibrating the: loading frame before taking photographs of the 
fringe pattern. Results of tests are given in tables III 
and IV and • on . f igur e s • 15 16 , and 17. 

Figure 15 indicates that experimental results on Model 
13 is in- good agreement with the theoretical value given by 
equation (22) for e = 0. This shows also that the choice 
of the original value of M oo; . does not affect the- plot of 

M or /M Q provided the values of M QC are large enough to 

produce sufficient fringes for accurate determination of 
nj - ng . 
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Figures 16 and 17 shov/ the comparison of experi- 
mental values of' M oc /M c for Models 14, 15 , and 16 with 

value's 'coinput ed from equation (24), Values of e on this 
curve were computed from measured values of y 1 and y £ , 
Table V gives values used for plots of solid lines on 
figures 16 end: 17. 

Table VI g i ve s ’ value s of P cr and y avv obtained by 

substitution of experimentally determined values of M 0( ,/ M Q 

in equation (25). These may be compared to .tjjQ theoretical 

tt I ' . 

values of P cr = — and the measured eccentricities on 

the models. I . 

Cjinolujs j.pn.s 

1. For beam columns with constant ; SI and with pin 
ends, the use of equation (25) for extrapolating experi- 
mental results appear to give good accuracy. 

" ’ ' ... . s . * * 

2. Range of M oc /M c values from 1.0 to appr oximat e ly 
0.25 should be covered by the tests, 

3. It is net necessary to have zero eccentricity of 
the end load. Results. can be analyzed by use of equation 
(25) to find the eccentricity present in the test. If 
values of M 0C /J4 C . are then desired f or other eccentricities, 

these may be obtained by use of equation (25) for any other 
arbitrary eccentricity. 

4. Eccentricities will have no effect on the critical 
load, but may have very rlo.rgd ef f es t - on the bending moments 
at values of axial load; les:s' -.than . the critical lead, If 

the bending, stress.es in the: b-a'ani i r e a c h t h'£ ‘-y i eld stress at 
values of. P less than .-the- critical load the allowable value 
of the axial load P wi-11 be very greatly affected by small 
variations in eccentricity. 

,5. On t he model the- lateral load’ should be sufficient, 
to give, values of M oc which are capable of accurate de— 

termination experimentally. This means that - n_ n-g 

should be at least 6. 
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TESTS OF CIRCULAR RING MODELS 
Purpo se 


The purpose of this series of tests was to give an ex- 
ample of the use of the basic similarity equations, (4) and 
(5/, and to determine whether continuous rings are struc- 
tures in which the a s sumpt i on of linearity of bending mo- 
ment with load is applicable. The bending moment, on a ring 
loaded with a diametrical .load P, is usually taken as 
(references): 

M = PR ( 0.3183 - 1 sin © ) (28) 


where 0 is the angle measured from a loaded diameter. At 

9 = ±90° = , this gives: 

2 



PR 


-0.1817 


This formula indicates no variation of moment with va- 
riation in the ratio R 4 /I or with variation in the ER 2 /P 
ratio. However, if the deflection of the ring has any ef- 
fect on the bending moment , .one or both of these ratios may 
be si giii f i cant . If the ratio ER 2 /P is significant, h 
variation in M^/g/PR would be' . expected with variations in 

that parameter. If. the moment of inertia of t*he ring is im- 
portant, a variation in M V ^ S /PR with R 4 /I 'would be ex- 
pected. For the circular ring as loaded here, other ratios 
as specified in equation (5) are satisfied with the excep- 
tion of p. , which is not significant in this problem. 


Models, Tests, and Computations 

Three circular ring models were constructed to dimen- 
sions shown in figure 18. These were cut from a sheet of 
polished and annealed Bakelite BT 61-893, 0.278 inch thick. 
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In table VII are listed tests run on the circular ring 
models. The average radius of the ring is denoted by R. 
Figure 19 shows a typical fringe photograph of a loaded ring. 

The test covered n range of ER e /P from 0.71x10 to 
12 . 6 Ox 10 5 and a range of R 4 / I from 2.99xl0 4 to 22.2X10 4 . 

Bending moments in circular rings may be computed with 
fair accuracy by using the same methods as discussed for 
straight beams. Since, however, the bending stress is not 
distributed lineerly in curved beams, the maximum bending 
moments "at the point 0 = ±90" were computed, for greater 
accuracy, as follows: 

1. The distance from the edge of the ringe to each 
fringe was measured on the fringe photograph. These measure- 
ments were then adjusted to model scale. Measurements were 
made f or both 0 = + 90° and —90°. Distances to fringes were 
averaged, pt the-s e two stations. The fringe order was plotted 
versus the average distance from the edge of the model, A 
typical plot for Model 17 under a 25— pound load is shown in 
figure 20. The values are recorded in table VIII. 

2, By planimeter integration of the area under the plot 
of n versus distance from' edge of model, the average fringe 
order for the beam was computed. The neutral axis was then 
the point at which the average fringe order was equal to the 
fringe order of the curve of n versus distance from edge 

of model/ 1 

3, At various distances y from the neutral axis, 
the fringe orders were measured. In table IX are shown 
these valuer for Model 1\7 under 25— pound lead. Values of 
y then were multiplied by fringe order n\. The product 

ny then, was plotted., Figure 21 shows the plot for" Model 17, 
load 25 Rounds. The area under, the curve of ny versus y 
then was measured. This is designated .as .5 . 

4. The bending iaotent across the. section under con— 1 
siderati bh may be written as, 


a 1 y b dy ' . ’ . (2 2) 



Mtr/e 
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where 


■' 



■ f 

’ \ % c 

** '*• : ; ' 



normal stress 

on section 



V. r . ' ■ ' . 


y 

distance 

f rou 

neutral 

axis 





b 

thicknes s 

of 

ring 






No 

..distance 

f r ora 

neutral 

axis 

t 0 

out br 

edgd of 

ring 

• ’ - 

yix T. 

i 

f.di st ance 

f. r oia 

neutral 

axi s 

t 0 

inner 

r : -;'x*T •:. * 

edge 7 of 

ring 


F#omI equat i 6fl (13), since a = 0,a x = -kn , equation^ 
( 29 ) r 'i?ec . qmes - r ... . > ;• 


^0 

n y dy 


« - w '■ ;• 
•v : 


M tr/ 2 = 


* - kb 


4 


i n i • i • o ! 
ni 

Therefore f 





k b S 



( 30 ) 


For mirror e rcj.etai led . di sens si on of this method of deter- 
mining tending moments, see reference 5. 


Re suit s 

• : In table- VII the values of M y ./PR are given for all 

• r . ~ T ' . ^ • h • ‘ 

tests. ' 'An ^examination of these results shows that for Model 
17* "the results are ;vqry nearly those given by equation (,28);. 
For Model's 'Iff and 19, however, the values of M y /PR are,, 

- ' * ' c5 

somewhat higher. In figure 22 are plotted average ya’^'asj . „Q.f 
M ' / -/PR versus . ER*/P for constant values of' the .ratio T - 

R /I. These curves show no significant variation with ER'/F # 

The variation’ that is shown is within experimental limits of 

accuracy. The curve of M , /PR versus H ^ / 1 , however, 

v 
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indicates an increase in the bending moment with an increase 
in R 4 /I, At low values of the ratio R 4 /l the value of 

/TB. checks almost exactly the value given by equation (28). 


Conclusions 

The photoelastic results show that the bending moments 

on circular rings are a function of the ratio R“/l. It is 
recommended that models of bulkhead rings and other similar 

structures be constructed so that the ratio., R 4 /l is the 

same on the model and full-scale structure. It also would 

be well to test the model at the same ratio of ER^/P on the 
model as for the full-scale structure if this is convenient. 

Even though little variation with ER C /P was shown in tests, 
it may be that at other values of the ratio, or for other 
types of loading, this would be significant. 

The photoelastic tests of these rings indicate that 
formulas in use for the calculation of bending moments in 
circular rings are in error for rings of large R 4 /I ratios. 
Furthermore these errors lead to nonconservative values. 

ANALYSIS ' OF STATICALLY INDETERMINATE FRAME 

, r "» , ^ « 

Purpose 


In statically indeterminate strictures, the loads and 
bending moments in various members of the. structure are func- 
tions of the relative stiffnesses of the members. If all 
members are loaded with axial loads that are small com- 
pared to their critical buckling loads, the analysis. of .the 
structure may be accomplished analytically by several 
commonly used methods. In structures where one or more of 
the members is loaded with an axial, load approaching its 
critical load, dif fifep.lt ies in the analysis occur due to the 
variation in the effective stiffness w h i c h ; t a k e's place in 
the axially loaded member as it approaches its critical load. 
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The photoelastic method, is useful in analyzing such a struc- 
ture, provided the model is constructed and loaded as speci- 
fied by the similarity principles previously discussed. 

The information desired for the stress analysis of the 
structures is as follows: r.-. 

1 # The total load at which buckling of the column mem- 
ber ..occurs 

2. The bending moment at every point in the structure 

3. The axial loads in all members of the structure 

4. The shear in the members 

If this information is obtained, unit stresses in the full- 
scale structure may be computed. 

If it were not for the effect of variation of -effective 
stiffness in the column member as its load is increased, the 
internal forces, reactions, bending moments, and so forth, 
would be linear functions of the external loads and these 
could be obtained for any convenient value of the applied 
load, say P x . The corresponding moments, internal forces, 
and so forth, in the various members could, then be found at 
any load P 2 by multiplying observed values at any load 
P x ..by the ratio, of P g /P 1 , 

When, however, one member of the structure is stressed 
by an axial load approaching its own critical load, the dis- 
tribution of the external loads among the members changes. 

The stresses in the individual members are then net linear 
functions of the external loads. In these structures the 
load on the model must be selected sc that the internal load 
distribution at a desired full-scale loading condition will 
be obtainable from the model results. In a structure of this 
type similarity conditions as specified in equation (8b) are 
applicable. 

Model and Tests 

To illustrate the use of the photoelastic method, two 
models were constructed as shown in figure 23 # 'These were 
constructed of polished and annealed Bakelite BT 61—893. 
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Tests were started on Model 20, tut the model failed 
before tests were completed. The results of these tests 
are not reported. Tests on Model 21 were run as follows: 

1. Fringe photographs of the central portion of the 
vertical leg were taken with the load P on the frame vary- 
ing from 40 to 100 pounds. Fringe photographs are shown in 
figure 24. 

2. Fringe photographs of the entire frame were taken 
at a load P = 40 pounds. These are shown on figure 25. 

3. Fringe photographs of the entire frame were taken 
at a load, of P = 75 pounds. These are shown on figure 26. 


Computations and Results 


The fringe patterns of figure 24 were analyzed in a 
manner .similar to that described for the beam— column models. 
Bending moments at a point 3.80 inches from -the pin point 
in the vertical member were determined. These are given in 
table X. .These moments correspond to the moment M c in 
e quat i on ( 25 ) . 

In order to analyze the results of this test, equation 
(20) is rewritten to conform to the particular loading con- 
dition on this. model. If the pin end of the vertical mem— 

- ber is considered to be equivalent to the right end of the 
span shewn .in figure 11, then, 


n = y fc = 6 « = b =0 



Equation (20) then becomes, 

, P 


M 


3 - P 


[K, - P'6,] | 


r l - 


a 


1 

2L 


( 31 ) 


P 1 is the axial load in the vertical member. 


where 
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','T , :*5 . 

Let 

P « 

= <t> x P 


L: • 


tj t }■, 

Mi 

:= 4J& 

t '> <1 
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' ‘ -M 


.* .-} J- 


.i 


Then, / s ince 5 : . , ; a } , * and <J> 2 , may 1 i be . f unc t i pn s . pf P, 

equation (31) may be written, ' ? £ J l 

* C 1 v ... " -* • ' - , • 5 


r * 


■7 r. ;:> 


. * x 


M c 



p *. 4> x p \ L / 


't v: . - f r: 


or 


/ h \ . 3 .. : 

P ( l a ) re P 

? — J— - _A f ( p ) 


Mr 


••A! & 


"*'(32) 


" ' •* 


n i: 


TH§ function .f(P) is not, in general, determinable 
analytically. An examination of the physics Cf this, how- 
ever, indicates that this function cannot he zero for'any 

' " ; A . ■ •; .1 : ' /. H. V ' - r ■ '<>•:•• 

= < . P \.*a “# I* : i h ' ^ 

value of ; P, -Consequent ly , if 

*• X i - t . w „ 1 

• •• '• . r. ■ u .- 1 ».• •. _ . 


•: *v- j: *x 

\ 




: f * 


against ^ P, this will -he zero w 


when 


is plotted 

M c ‘ ' . 7 .-: V-,.,,. 

n " I 6,. 0 

- — — P , is* equal to 

<P . ’ f *t . • j . > 


zero, *Kie value of P. .at this condition is then the' value 
of the total load on the structure which will tend pro- 
duce an infinite bending moment in the vertical l*egy- : ‘ Tjhi s 
may be* ^'called the cribircal value of .P for t be' vertical 
column. The critical value of P for a strucfur e of tliis 
type m'ay be interpreted a.s the total load on the structure 
at which the buckled member will carry no additional axial 
load even if P is further increased. 

\ ‘ i J l '• 

1 ^ 

Z 2 J 

In figure 27. 1 is plotted against P. 

M c 

The curve through experimental points was extrapolated to 
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/ 

the point 

total load of 122 pounds, the vertical member is essenti- 
ally buckled. Since the structure is redundant, this does 
not necessarily mean that this is the ultimate load on the 
structure. It indicates, however, that any additional load 
will have to be carried entirely on the left support of the 
horizontal member. Consequently, at loads above 122 pounds 
the -bending moments in the left portion of the •beam will in- 
crease much more rapidly than at loads less than 122 pounds. 

To obtain the bending moment at every section of the 
structure at given loads full— s-cele, the e qu i’ vr. 1 % rot loads 

are computed so that P/EL*" for the model and the full- 
scale structure are identical. For the tests of Model 21, 
assume that these model loads have been determined to be 
40 and 75 pounds. 

The first step in the analysis was to determine the 
fringe orders at various sections along the boundary of the 
model. These were identified by observing the formation of 
the fringes as the load w as increased. These fringes were, 
recorded on figures 25 and 26 apd in. tables XI and Xll, At 
sections of; maximum moment, it impossible to plot fringe 
order against depth of beam and ' to*use faired values’ of* 
fringe orders at the boundary. 'At sections where bending 
moments are not maximum,' this method is not applicable. 

Good results can be obtained, however, by using " : the esti- 
mated values- ‘of f r i nge ~ot d er s at the outer boundaries and 
then determining the bending moments from equation (14). 

This was done at various sections of the horizontal and 
vertical members of Model 21 for loads of both 40 and 75 
pounds. These bending moments are plotted also’ on figures 
25 and 26. -The shear forces in the members may be obtained 
by measurement of the slopes of the moment curves. 

The axial loads in the members may be obtained by 
taking the average fringe order across the section. The 
axial load is then given by J 


h, 


M r 


= 0. This indicates that, at a 


p 

axial 




n av th 
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where t is the thickness 'and- ~ h the width of the cross 
section. The accuracy of determination of n aV is poor 

if the average fringe arder i s: 7 s.ma ; .l 1 . Usually, however, a 

relatively large error in this axial load does net affect 

the combined stress, . . - >. ’ ...” . 

* .. • «.» 

, Mh 

a = 

to any great extent. If a is large, then it may be 
determined with greater accuracy since the average fringe 
order will be larger. . • : ... . 

In the model tested, the axial. load in the column is 
due to shear carried in the beam. At the 40-pound external 
load the axial load in the column obtained by considering 
the shear in the right-hand portion of the beam was found 
to be approximately 18.0 pounds. At the 75— pound external 
load the load carried tiy the column was found tp be approxi- 
mately 34,0 pounds. It is seen that the four .desired items 
of information for the full-scale structure are obtainable. 

In beam— column structures, part icularly .when the beam- 
column member is redundant, it is very difficult, by usual 
methods, to determine allowable loads.- In structures with 
linear behavior, if a factor of safety on the ultimate unit 
stress of 2.0 is desired, it is necessary merely to determine 
the external load giving a unit stress equal to one— half of 
the ultimate stress for the material and to consider this 
the design load. In beam— column structures, however, ex- 
ternal loads which produce bending stresses equal to one- 
half the allowable stresses actually may be loads very near 
the external load at which the structure will fail by buck- 
ling. Also as axial- load's in beam— column members approach 
critical values, the bending moments increase much more 
rapidly than the load. Hence ultimate stresses may be 
reached at external Loads very little greater than those 
producing only one-half the ultimate stress at the most highly 
stressed section of the beam column. In the case of long 
columns, of course, loads just slightly less than the criti- 
cal may produce stresses only a small fraction of the ulti- 
mate stress. Tests of photoelastic models may be used to 
indicate the stress conditions at various loads, and an accu- 
rate determination of the allowable loads for the full-scale 
structure may be made. 
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CONCLUSIONS 


By intelligent application of the similarity princi- 
ples given, photpelastic models nay be ccnstruqtsd to repre- 
sent many types of structures that cannot be... r ead i ly analyzed 
by analytical methods. These include: 

1. Beam columns of varying El 

2. Structures in which def ormat i ons are large enough 

sc that their effect on stress distribution is 
not negligible (such as flexible circular rings) 

3. Statically indeterminate structures in which in- 

ternal forces, moments, and so forth, are, not 
linear functions of the external loads 


T'he examples given indicate that the techniques of test 
and meth’ads of analysis used give accurate values, o.f critical 
loads, bending moments, and so forth, for the p h.o/t cel a s t i c 
models. Corresponding full-scale values may ; be accurately 
obtained if similarity principles have been ^applied properly. 


Oregon State College, 

Corvallis, Ore., June 1945. 
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TABLE I 

Tests and Results on Models No. 1 to 11 


Photo 

No. 

Model 

No. 

— 

Test 

No. 

P 

Q 

“b^t 

M 

oc 

M c 

Computed 0 

P cr 

1-B 

1 

1 

0 

26.15 

8.56 

1.000 

1200 

2 

1 

1 

0 

52.30 

15.90 

1.000 

1200 

3 

1 

1 

75 

52.30 

17.88 

.890 

1200 

4 

1 

1 

150 

52.30 

19.60 

.812 

1200 

5 

1 

1 

225 

52.30 

21.38 

.744 

1200 

6 

1 

1 

275 

52.30 

22.50 

.707 

1200 

7 

1 

1 

300 

52.30 

23.20 

.685 

1200 

8 

2 

1 

0 

26.15 

10.86 

1.000 

835 

9 

2 

1 

75 

26.15 

12.70 

.855 

835 

10 

2 

1 

150 

26.15 

14.52 

.748 

835 

11 

2 

1 

225 

26.15 

16.70 

.651 

835 

12 

2 

1 

200 

26.15 

16.10 

.674 

835 

13 

2 

1 

225 

26.15 

16.70 

.651 

835 

14 

3 

1 

0 

20.92 

11.70 

1.000 

530 

15 

3 

1 

75 

20.92 

14.58 

.803 

530 

16 

3 

1 

150 

20.92 

16.60 

.785 

530 

17 

3 

1 

225 

20.92 

19.60 

.596 

530 

18 

3 

1 

275 

20.92 

23.38 

.501 

530 

19 

3 

1 

300 

20.92 

26.50 

.442 

530 

20 

4 

1 

0 

12.55 

9.90 

1.000 

294 

21 

4 

1 

50 

12.55 

13.25 

.747 

294 

22 

4 

1 

100 

12.55 

16.34 

.606 

294 

23 

4 

1 

125 

12.55 

18.50 

.535 

294 

24 

4 

1 

150 

12.55 

21.30 

.465 

294 

25 

4 

1 

175 

12.55 

25.48 

.389 

294 

26 

5 

1 

0 

7.32 

9.30 

1.000 

146 

27 

5 

1 

25 

7.32 

11.34 

.820 

146 

28 

5 

1 

50 

7.32 

14.10 

.659 

146 

29 

5 

1 

75 

7.32 

19.02 

.488 

146 

30 

5 

1 

90 

7.32 

24.20 

.384 

146 

31 

5 

1 

100 

7.32 

29.28 

.317 

146 

32 

6 

1 

0 

2.09 

5.80 

1.000 

45.7 

33 

6 

1 

10 

2.09 

7.00 

.828 

45.7 

34 

6 

1 

15 

2.09 

8.20 

.707 

45.7 

35 

6 

1 

20 * 

2.09 

9.80 

.592 

45.7 

36 

6 

1 

25 

2.09 

11.30 

.513 

45.7 

37 

6 

1 

30 

2.09 

15.20 

.381 

45.7 


°By interpolation of tables in reference (2) 
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TABLE I (Cont’d) 

Tests and Results on Models No. 1 to 11 


Photo 

No. 

Model 

No. 

Test 

No. 

P 

Q 

n B~ n T 

“ E | f 

r^rl 

Jompu'ted 

*cr 

38 

6 

1 

0 

4.18 

11.30 

1.000 

45*7 

39 

6 

1 

10 

4.18 

14.28 

.792 

45.7 

40 

6 

1 

20 

4.18 

18.50 

.612 

45.7 

41 

6 

1 

25 

4.18 

25.80 

.438 

45.7 

42 

6 

1 

30 

4.18 

32.90 

.343 

45.7 

43 

7 

2 

0 

26.15 

10.24 

1.000 

401 

44 

7 

2 

75 

26,15 

12.30 

.832 

401 

45 

7 

2 

150 

26.15 

16.00 

.640 

401 

46 

7 

2 

225 

26.15 

19.59 

.523 

401 

47 

7 

2 

275 

26.15 

23.00 

.445 

401 

48 

7 

2 

300 

26.15 

25.58 

.401 

401 

49 

8 

2 

0 

15.70 

8.88 

1.000 

309 

50 

8 

2 

100 

15.70 

12.34 

.720 

309 

51 

8 

2 

175 

15.70 

16.16 

.550 

3C9 

52 

8 

2 

200 

15.70 

18.84 

.472 

309 

53 

8 

2 

225 

15.70 

23.56 

.377 

309 

54 

8 

2 

250 

15.70 

33.20 

.268 

309 

55 

9 

2 

0 

12.55 

9.54 

1.000 

216 

56 

9 

2 

50 

12.55 

12.68 

.752 

216 

57 

9 

2 

100 

12.55 

16.70 

.570 

216 

58 

9 

2 

125 

12.55 

20.38 

.468 

216 

59 

9 

2 

150 

12.55 

26.30 

.362 

216 

60 

9 

2 

175 

12.55 

41.00 

.232 

216 

61 

10 

2 

0 

8.37 

10.65 

1.000 

123 

62 

10 

2 

15 

8.37 

12.25 

.869 

123 

63 

10 

2 

25 

8.37 

13.05 

.815 

123 

64 

10 

2 

50 

8.37 

16.15 

.659 

123 

65 

10 

2 

75 

8.37 

22.45 

.474 

123 

66 

10 

2 

100 

8.37 

44.40 

.239 

123 

67 

11 

2 

0 

4.18 

10.50 

1.000 

45.7 

68 

11 

2 

10 

4.18 

14.00 

.750 

45.7 

69 

11 

2 

20 

4.18 

18.40 

.572 

45.7 

70 

11 

2 

25 

4.18 

24.35 

.432 

45.7 

71 

11 

2 

30 

4.18 

31.50 

.333 

45.7 

72 

11 

2 

35 

4.18 

43.25 

.243 

45.7 


°By interpolation of tables in reference (2). 


TABI£ II 

Analysis of Eccentricities and Critical Load 


Model 

p i 

1 

P 2 

2P 1 P 2 

3P 1 

3 V P 2 

« 

/S=P 
~ cr 

Q 

M oc 

M oc 
•» ■ ■■■— 

p i 

^oc 

P or 

+ 

^av 

% error 
P cr 

1 

475 

805 

765000 

1425 

620 

1232 

52.30 

78.50 

- 0.1651 

0.1102 

- 0.0549 

2.9 

2 

330 

563 

371000 

990 

427 

875 

26.15 

39.20 

- .1187 

.0895 

- .0292 

4.8 

3 

265 

400 

212000 

795 

395 

536 

20.92 

31.40 

- .1182 

.1171 

- .0011 

1.2 

4 

134 

213 

57100 

402 

189 

302 

12.55 

18.80 

- .1402 

.1245 

- .0157 

2.8 

5 

73 

109 

15920 

219 

110 

145 

7.32 

11.00 

- .1503 

.1518 

.oon 

- .9 

6 

25 

36.2 

1810 

75 

38.8 

46.6 

2.09 

3.14 

- .1255 

.1348 

.0093 

2.0 

6 

23 

34.8 

1602 

69 

34.2 

46.8 

4.18 

6.27 

- .2720 

.2680 

.0040 

2.4 

7 

252 

338 

170200 

756 

418 

407 

26.15 

39.20 

- .1558 

.1922 

.0434 

1.6 

8 

190 

257 

97700 

570 

313 

312 

15.70 

23.60 

- .1240 

.1510 

.0270 

1.0 

9 

3.20 

170 

40800 

360 

190 

214 

12.55 

18.82 

- .1568 

.1760 

.0192 

- .9 

10 

70 

99.0 

13380 

210 

111 

125 

8.37 

12.55 

- .1792 

.2008 

.0216 

1.7 

11 

23.0 

34.5 

1585 

69 

34.5 

45.8 

4.18 

6.27 

- .2730 

.2740 

.0010 

.2 


°Fr om Eq uati on 26 
•♦From Equation 27 


TABLE m 

Results of Tests on Model 13 


Photo 

Model 

Test 

P 

Q 

V n T 

"oc 

Computed 0 

No. 

No. 

No. 



M c 

P cr 

73 

13 

3 

0 

3.14 

5.50 

1.000 

98.2 

74 

13 

3 

30 

3.14 

7.92 

.695 

98.2 

75 

13 

3 

50 

3.14 

11.10 

.496 

98.2 

76 

13 

3 

60 

3.14 

13.78 

.399 

98.2 

77 

13 

3 

70 

3.14 

18.06 

.305 

98.2 

78 

13 

3 

75 

3.14 

22.12 

.249 

98.2 

79 

13 

3 

0 

6.28 

11.00 

1.000 

98.2 

80 

13 

3 

30 

6.28 

15.75 

.698 

98.2 

81 

13 

3 

50 

6.28 

22.42 

.490 

98.2 

82 

13 

3 

60 

6.28 

29.00 

.379 

98.2 

83 

13 

3 

70 

6.28 

38.50 

.286 

98.2 

84 

13 

3 

75 

6.28 

45.50 

.242 

98.2 

85 

13 

3 

0 

9.41 

16.50 

1.000 

98.2 

86 

13 

3 

30 

9.41 

23.75 

.694 

98.2 

87 

13 

3 

50 

9.41 

33.80 

.488 

98.2 

88 

13 

3 

60 

9.41 

42.50 

.388 

98.2 

89 

13 

3 

70 

9.41 

54.10 

.305 

98.2 


op- 

P cr- 12 
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TAB1E IV 

Tests and Results on Models 14, 15, 16 


Photo 

Model 

P 

Q 

V”t 

M 

oc 

i 

H 





M c 

cr calc 




.100", e = 

- 3.21 



90 

14 

0 

15.69 

9.29 

1.000 

0 

91 

14 

125 

15.69 

16.95 

.548 

.167 

92 

14 

225 

15.69 

26.00 

.357 

.301 

93 

14 

325 

15.69 

38.95 

.239 

.435 



y aT =< 

).100", e = 

: 3.21 



94 

14 

0 

15.69 

8.40 

1.000 

0 

95 

14 

125 

15.69 

5.40 

1.556 

.167 

96 

14 

225 

15.69 

1.70 

4.941 

.301 

97 

14 

450 

15.69 

-15.20 

- .553 

.602 



- °' 

.046", e = 

- 1.29 



98 

15 

0 

10.46 

8.85 

1.000 

0 

99 

15 

100 

10.46 

14.55 

.608 

.225 

100 

15 

175 

10.46 

21.70 

.408 

.394 

101 

15 

275 

10.46 

38.80 

.228 

.618 



y = ( 

3.046", e = 

: 1.29 



102 

15 

0 

10.46 

8.57 

1.000 

0 

103 

15 

100 

10.46 

7.78 

1.102 

.225 

104 

15 

200 

10.46 

7.46 

1.148 

.451 

105 

15 

300 

10.46 

5.95 

1.440 

.675 



^av O ' 

,005", e = 

-0.11 



106 

16 

0 

7.32 

8.97 

1.000 

0 

107 

16 

60 

7.32 

12.04 

.745 

.251 

108 

16 

100 

7.32 

15.42 

.582 

.418 

109 

16 

160 

7.32 

27.53 

..326 

.669 



^av - ( 

3.005", e = 

= 0.11 



110 

16 

0 

7.32 

9.28 

1.000 

0 

111 

16 

75 

7.32 

14.05 

.660 

.313 

112 

16 

125 

7.32 

18.92 

.490 

.522 

113 

16 

175 

7.32 

32.08 

'.289 

.732 


°P calc = 
cr 


TT^EI 

ft 2 
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.. , “oc 1 -A 

Values from: r , — = z — 

M 1 -jie 


<4 

M /M 
oc' c 

e = -3.21 

M„ /M 
oc' c 

e = 3.21 

M /M 
oc' c 

e = -1.29 

M Al 
oc' c 

e = 1.29 

P 

< 0( ? 

a> 

M /M 
oc c 

e = 0.11 

.00 

1.000 

1.00 

1.000 

1.000 

1.000 

1.000 

#05 

.820 

1.13 

.892 

1.015 

.945 

.955 

.10 

.682 

1.32 

.797 

1.032 

.889 

.910 

.15 

.575 

1.63 

.712 

1.053 

.836 

.864 

.20 

.488 

2.22 

.637 

1.078 

.782 

.818 

.25 

.416 

3.75 

.565 

1.123 

.730 

. 772 

.30 

.357 

17.50 

.505 

1.142 

.678 

.724 

.35 

.307 

- 5.41 

.447 

1.185 

.627 

.675 

.40 

.263 

- 2.14 

.396 

1 . 2 U 

.574 

.628 

.45 

.225 

- 1.25 

.347 

1.312 

.523 

.578 

.50 

.192 

- .83 

.304 

1.408 

.473 

.529 

.55 

.163 

- ..59 

.263 

1. 543 

.423 

.479 

.60 

.137 

- .43 

.225 

1.771 

.375 

.428 

.65 

.114 

- .32 

.191 

2.160 

.326 

.377 

.70 

.093 

- .24 

.157 

3.090 

.279 

.325 

.75 

.073 

- .18 

.127 

7.590 

.231 

.273 

.80 

.056 

- .13 

.098 

- 6.140 

.183 

.219 

.85 

.040 

- .09 

.072 

-1.550 

.137 

.165 

.90 

.026 

- .05 

.046 

- .620 

.091 

.111 

.95 

.016 

- .02 

.023 

- .220 

.045 

.056 

1.00 

.000 

.00 

.000 

.000 

.000 

.000 


TABLE VI 

Values of P and y Obtained from Tests 
cr av 

on Models 14* 15* and 16 


Model 

No. 

Measured 

Eccentricity 

Experimental 

M /M 
oc 7 c 

Load P 
lb 

°P 

cr 

lb 

°y 

y av 

in. 

P cr 

Theoretical 


y inches 




14 

- 0.100 

0.357 

225 

770 

- 0.102 

748 

14 

- 0.100 

.239 

325 




14 

0.100 

4.9a 

225 

745 

0.090 

748 

14 

0.100 

- .553 ■ 

450 

15 

- 0.046 

0.408 

175 

480 

0.050 

445 

15 

- 0.046 

.228 

275 

15 

0.046 

1.148 

200 

468 

0.040 

445 

15 

0.046 

1.440 

300 

16 

- 0.005 

0.582 

100 

236 

0.001 

239 

16 

- 0.005 

.326 

160 

16 

16 

0.005 

0.005 

0.490 

.289 

125 

175 

247 

- 0.001 

239 


•Obtained by substitution of load P and experimental M qc /M c in Equation 25# 
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TABLE VII 

Circular Ring Tests 
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Fhoto 

No. 

Model 

No. 

Test 

No. 

Load 

R 

PR 1 

Va 

M rr/2 

PR 

ER 2 

P 

R 4 

I 

114 

17 

7 

12.5 

2.30 

28.8 

- 5.24 

-0.182 

2.82* 10 5 

2.99*10 4 

115 

17 

7 

25.0 

2.30 

57.5 

-10.62 

4wm 

.184 

1.41* 10 5 

2.99 X 10 4 

116 

17 

7 

37.5 

2.30 

86.3 

-15.79 

- 

.183 

,94*10 5 

2.99 X10 4 

117 

17 

7 

50.0 

2.30 

115.0 

-20.58 

- 

.180 

.71X10 5 

2.99 X10 4 

118 

18 

7 

5.0 

2.35 

11.8 

- 2.33 

» 

.198 

7.36X 10 5 

7.93* 10 4 

119 

18 

7 

10.0 

2.35 

23.5 

- 4.35 

- 

.185 

3.68X10 5 

7.93 X10 4 

120 

18 

7 

13.0 

2.35 

30.6 

- 5.86 

- 

.192 

2.84X10 5 

7.93 X 10 4 

121 

18 

7 

16.0 

2.35 

37.6 

- 7.21 

•a 

.192 

2.30 *10 5 

7.93 *10 4 

122 

18 

7 

18.0 

2.35 

42.3 

- 7.90 

- 

.187 

2.05X10 5 

7.93 *10 4 

123 

18 

7 

20.0 

2.35 

47.0 

- 8.84 

- 

.188 

1.84X10 5 

7.93 X 10 4 

124 

18 

7 

23.0 

2.35 

54.0 

-10.10 

9WC 

.187 

1.60X.10 5 

7.93 X 10 4 

125 

18 

7 

25.0 

2.35 

53.7 

-10.91 

mm 

.186 

1.47X10 5 

7.93 *10 4 

126 

19 

7 

3.0 

2.38 

7.15 

- 1.485 


.208 

12.60X10 5 

22.20X10 4 

127 

19 

7 

6.0 

2.38 

14.30 

- 2.880 

- 

.202 

6.30X10 5 

22.20*10 4 

128 

19 

7 

9.0 

2.38 

21.40 

- 4.350 

- 

.203 

4. 20 X10 5 

22.20* 10 4 

129 

19 

7 

10.5 

2.38 

25.00 

- 4.974 

- 

.199 

3.60X10 5 

22.20* 10 4 

130 

19 

7 

12.0 

2.38 

28.55 

- 5.700 

- 

.200 

3.15X10 5 

22.20 X10 4 


TABLE VIII 

Fringe Measurements for Model 17, P s 25*0 lb 


Fringe order 
n 

Distance from outer 
edge of photograph 
in inches 

Distance from outer 
edge of model. (Model 
scale), inches 


0 = - 90° 

9 — 90® 

Average 


4 

0.035 

0.025 

0.030 

0.041 

3 

.053 

.055 

.054 

.073 

2 

.078 

.083 

.080 

.108 

1 

.104 

.107 

.105 

.142 

0 

.127 

.132 

.129 

.174 

- 1 

.149 

.155 

.152 

.205 

- 2 

.170 

.175 

.173 

.232 

- 3 

.193 

.197 

.195 

.263 

- 4 

.213 

.216 

.214 

.289 

- 5 

.231 

.234 

.232 

.313 

- 6 

.248 

.253 

.250 

.338 

Inner Edge 

.259 

' .259 

.259 

.350 
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TABLE IX 

Calculations for Model 17, Load = 25.0 lb 
(Values of y and n from Figure 20) 


Distance from 
neutral axis 

y 

Fringe Order 
n 

ny 

0.183 

5.41 

0.990 

.143 

4.30 

.615 

.103 

3.16 

.326 

.063 

2.08 

.131 

.023 

.75 

.017 

.000 

.00 

.000 

- .037 

- 1.25 

.046 

- .077 

- 2.70 

.208 

- .117 

- 4.15 

.485 

- .157 

- 5.78 

.908 

- .167 

- 6.20 

1.035 


TABLE X 

Bending Moments in Section 3.80 in. from Pin Joint 
in Vertical Member of Model 21 


Photo 

No. 

Model 

No. 

Test 

No. 

Load 

P 

M c 

h- a 

p<V 1) 

Mt 

151 

21 

8 

40 

0.875 

347 

152 

21 

8 

50 

1.305 

291 

153 

21 

8 

60 

1.750 

260 

154 

21 

8 

65 

2.136 

231 

155 

21 

8 

70 

2.410 

220 

156 

21 

8 

75 

2.940 

194 

157 

21 

8 

80 

3.580 

170 

158 

21 

8 

85 

4.240 

152 

159 

21 

8 

90 

5.440 

126 

160 

21 

8 

95 

6.660 

108 

161 

21 

8 

100 

8.200 

93 


°P, ^ defined in Figure 23, 
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TABLE H 

Bending Moments on Model 21, P =40*0 lb 


Horizontal Beam 

Section 

Distance from 

n T" “B 

Bending 


left support 

1 D 

moment 

a 

0.00 



b 

.30 

11.2 

- 20.80 

c 

.62 

7.5 

- 13.95 

d 

.93 

3.5 

- 6.50 

e 

1.26 

0.0 

0.00 

f 

1.56 

- 3.1 

5.76 

g 

1.87 

- 7.0 

13.00 

h 

2.20 

- 10.4 

19.30 

i 

2.50 

- 13.6 

25.30 

3 

2.82 

- 16.9 

31.40 

k 

3.13 

- 20.0 

37.20 

1 

3.46 

- 16.0 

29.80 

m 

3.75 

- 13.0 

24.20 

n 

4.07 

- 9.9 

18.40 

0 

4.38 

- 6.6 

12.30 

P 

4.69 

- 3.6 

6.70 

q H 

5.28 




Vertical Column 

Section 

Distance from 

n T“ “b 

Bending 


pin joint 

i D 

moment 


7.87 


, 

r V 

6.55 

0.2 

- 0.069 

s 

5.91 

1.1 

- .377 

t 

5.30 

2.0 

— .685 

u 

4.65 

2.2 

- .755 

Y 

4.05 

1.9 

- .651 

TT 

3.41 

1.6 

- .549 

X 

2.79 

1.3 

- .446 

7 

2.15 

1.0 

- .343 

z 

1.54 

.9 

- ..308 

aa 

.90 

.6 

- .206 

bb 

.28 

— 

— — 
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TABLE XU 

Bending Moment on Model 21, P = 75.0 lb 



Horizontal Beam 


Section 

Distance from 


Bending 


left support 

1 D 

moment 

a 

0.00 

— — ... 


b 

.31 

23.2 

- 43.20 

c 

.63 

16.0 

- 29.80 

d 

.93 

8.5 

- 15.80 

e 

1.26 

1.0 

- 1.86 

f 

1.56 

- 5.2 

9.63 

g 

1.88 

- 12.5 

23.20 

h 

2.19 

- 19.0 

35.40 

i 

2.51 

- 25.3 

47.00 

i 

2.82 

- 33.7 

62.60 

k 

3.14 

— — 

— — 

je 

3*44 

- 33.3 

61.90 

m 

3.75 

- 27.7 

51.50 

n 

4.06 

- 22.1 

41.20 

0 

4.38 

- 15.9 

29.60 

P 

4.70 

- 10.0 
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Figure l.“ Dimensions of models I, 2, 3, 4, 5, 6 and 13. 
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Figure 2 - Dimensions of models 7,8,9 ,10 and II. 
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Figure 3.- Photograph of model 3 set up for test. Figure 5- Fringe photographs 

for model 6; P = 0 
to 30 lb, Q = 4.18 lb. 
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Figs. 4 



13 to 16. 



12 . 


X 




<r 


Figure 7.- Mohr's circle of stress at section of maximum moment 
in beam-column. 
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Figs. 8,9 



Figure 8 - Plots of M QC /M c against P for models 1,2,3 and 4. 



Figure 9 - Plots of M QC /M C against P for models 5,6 and II. 
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Figs. 12,13,14 
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Figure 12.- Dimensions of model 14. 
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Figure 13 - Dimensions of model 15. 
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Figure 14 - Dimensions of model 16. 
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Figs. 15,16 
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Figure 17 - Plots of M QC /M c against ct for models 14, 
15 and 16 with positive eccentricities. 
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Figure 18.- Dimensions of circular ring models 17, 18 
and 19. 
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Figure 19.- 


Fringe photograph for model 17; P 


25 lb. 






DISTANCE FROM OUTER EDGE OF MODEL, INCHES 

Figure 20.- Plot of fringe order against distance from outer edge 
of model at 0 = ±9O° for model 17, P = 25 lb. 
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Figure 21.- Plot of ny against y for model 17, P = 25 lb. 
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Figure 22.- Plots of NU/o/PR f° r circular ring models against 
ER 2 /P and R 4 /I. 
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Figure 23.- Dimensions of models 20 and 21. 
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Fig. 24 



Figure 24.- Fringe photographs of mid-section of 
vertical member of model 21. 
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BENDING MOMENTS POSITIVE IF UPPER EDGE 
OF BEAM OR RIGHT EDGE OF COLUMN IS IN 
COMPRESSION. 


NUMBERS AT SIDES OF FRINGE PHOTOGRAPHS 
GIVE FRINGE ORDERS AT SECTIONS MARKED. 
FRINGE ORDERS ARE POSITIVE FOR TENSILE 
STRESS AT BOUNDARY. 


MOMENTS IN BEAM, M x = -l.86Cn T -n^) 
MOMENTS IN COLUMN, My = -0.343(n T - n B )- 
WHERE'- 

n T = FRINGE ORDER AT TOP OF BEAM AND 
RIGHT SIDE OF COLUMN. 
n B = FRINGE ORDER AT BOTTOM OF BEAM 
AND LEFT SIDE OF COLUMN. 
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Figure 25.- Fringe pattern and bending moment diagrams for model 21, 
P = 40 lb. 
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BENDING MOMENTS POSITIVE IF UPPER EDGE 
OF BEAM OR RIGHT EDGE OF COLUMN IS IN 
COMPRESSION. 


NUMBERS AT SIDES OF FRINGE PHOTOGRAPHS 
GIVE FRINGE ORDERS AT SECTIONS MARKED. 
FRINGE ORDERS ARE POSITIVE FOR TENSILE 
STRESS AT BOUNDARY. 


MOMENTS IN BEAM, M x = -1.86 (ri T - n B ). 

MOMENTS IN COLUMN, Mw= -0.343 (n T -n 0 ). 
WHERE: 

n T = FRINGE ORDER AT TOP OF BEAM AND 
RIGHT SIDE OF COLUMN. 
n B = FRINGE ORDER AT BOTTOM OF BEAM 
AND LEFT SIDE OF COLUMN. 


— - 1.8 r 


0.6 S 


— 05 t 


1.2 U 


— 1.6 V 


— 1.6 W 


— 1.5 X 


— o.i y 


0.7 2 


— 1.5 aa 


bb 


Figure 26.- Fringe pattern and bending moment diagrams for model 21 
P= 75 lb. 
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Figure 27.- Variation in bending moment, in vertical member of model 21, with load. 
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